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Abstract. The energy input into the interstellar medium in Ultraluminous 
Infrared Galaxies (ULIRGs) is enormous, regardless of the nature of the 
power source. I discuss some of the major consequences for the structure 
and energetics of the ISM in these galaxies. Observationally, the column 
t — . densities in the nuclear regions of ULIRGs are known to be very high, 

which makes distinguishing starbursts from AGN quite difficult. The level 
of energy and momentum injection means that the pressure in the ISM must 
Q\ \ be extremely high, at least 3 — 4 orders of magnitude larger than in the local 

ISM or typical giant molecular clouds. It also means that the luminosity 
of GMCs in ULIRGs must be very high, as they must radiate many times 
their binding energy over their lifetimes. I briefly review the influence which 
X-ray irradiation can have on the ISM in AGN-powered ULIRGs. Finally, I 
show that the presence of PAH features in ULIRGs does not imply that they 
must be starburst-dominated, since at the column densities and pressures 
typical of the ISM in ULIRGs PAHs can survive even at tens of parsec 
distances from the AGN. 

1. Introduction 

The observed characteristics of ULIRGs (luminosities Ljr > 10 12 Lq, size 
scales r ;$ a few hundred pc) imply typical energy densities u ~ 10~ 8 erg 
cm~ 3 . This is 3-4 orders of magnitude larger than the typical energy den- 
sities (in gas, photons, and cosmic rays) in the interstellar medium (ISM) 
in the solar neighborhood (e.g., Jura 1987). In this review, I will discuss 
how both the production and release of such energies is expected to impact 
the ISM. I will concentrate on those aspects most relevant to ULIRGs, i.e., 
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very powerful starbursts and luminous active galactic nuclei (AGN). Due 
to space limitations, the discussion is necessarily terse. In particular, the 
referencing is by no means complete, and should be regarded largely as a 
pointer to the relevant literature (a good place to start is the recent review 
by Sanders & Mirabel 1996); my apologies to those whose contributions 
have of necessity been left out. 

2. Basic Parameters 

2.1. STARBURST: 

Bolometric Luminosity The frequency-integrated energy output from a 
starburst can be approximated as 

Lbo1 ~ 10l ° (i^M ' (rSo) (iSe) Lq (1) 

for a constant star formation rate (SFR) burst of age At, upper and lower 
mass cutoffs m u and mi, and a = 0.23 for mi < 1 Mq (Scoville & Soifer 
1991; see also Smith, Lonsdale, & Lonsdale 1998). (Here and frequently 
throughout this review, adopted scalings with mass, density, etc. should be 
regarded with the caution usually reserved for possibly rabid housepets.) 
This implies that powering ULIRGs by star formation requires 

modulo the details of the assumed initial mass function. Thus ULIRGs must 
have SFR in excess of 100 Mq yr _1 , barring extreme assumptions about 
the IMF. 

Lyman Continuum Photons The number of Lyman continuum (E > 13.6 
eV) photons produced per second by a starburst is approximately (Scoville 
& Soifer 1991) 

((3 = 2, m u < 45 M ; (3 = 1.35, m u > 45 Mq), where I have again assumed 
a constant SFR burst with an age of 10 8 years. Thus 

for a typical ULIRG. Since the radiation is stellar in origin, the spectrum 
of radiation is sharply cutoff above a few tens of eV. 
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Supernova Rate The implied supernova rate for a starburst-powered ULIRG 
is 

This rate is not very sensitive to the assumptions made about the IMF, 
since the predecessors of most SNe have M ~ 10 M & . 

2.2. AGN: 

Bolometric Luminosity The accretion luminosity for an accreting black hole 
(presumed, for many reasons, to be the power source in AGN) which is 
converting rest-mass energy to radiation with an efficiency e is 

L m = 1.5 x 10" (£) (-^) L (6) 

independent of the mass of the hole. However, the Eddington limit (the 
luminosity above which radiation pressure will halt accretion for steady, 
spherically-symmetric flow) is 

iEDD ,3.4 X 10»(^)(^) is . (7) 

If we require L acc < L EDD , we must have 

M,„ > 4 x HT (£) (-^j M S . (8) 

Lyman Continuum Photons The typical spectrum of an AGN is a power- 
law in energy, with a flux density oc v~ a , with a ~ 1 {i.e., the spectrum 
is flat, with approximately equal energy per decade of photon energy) . For 
unobscured (Type 1) AGN, the X-ray luminosity L x ~ O.lLbol, and the 
spectral index in the hard X-ray regime a ~ 0.7. The fraction of energy 
emitted in the 2-10 keV band decreases slowly with increasing luminosity 
(see the summary in Mushotzky, Done, & Pounds 1993). For a spectrum 
with vf u = constant from 1 eV to 100 keV, the number of Lyman continuum 
photons is 

essentially identical to the number produced by a starburst. Although an 
AGN spectrum extends to far higher energies than that of a starburst, most 
of the ionizing photons by number are between 13.6 and 100 eV, just as in 
a starburst. 
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3. Fueling 



3.1. STARBURST: 



Forming stars at a rate M s M & yr requires gas utilization at a rate 




(10) 



where e* is the efficiency of star formation, i.e., the fraction of the mass of 
a star-forming cloud which is actually transformed into stars in the cloud's 
lifetime. The observed value e* ~ 0.01 — 0.1 in Galactic molecular clouds. 
(Note, however, that in order to form a bound stellar system, such as a 
globular cluster, simple binding energy arguments show that the efficiency 
must be much higher, e* £ 0.4; if the so-called "proto-globular clusters" 
seen, for example, in NGC 4038/4039 [Whitmore & Schweizer 1995] are in 
fact bound systems, then the SFE must have been quite high.) Of this gas 
used by star formation, an amount 



is lost from the system. In this equation f^ s is the fraction of the non-stellar 
remainder of the star-forming cloud which has either been blown entirely 
out of the system (e.g., in a galactic superwind: see §4), or is either too hot 
or too dispersed to reform into star-forming clouds on the characteristic 
timescale of the burst; / rem is the fraction of the stellar mass that is in the 
form of stellar remnants or stars which have masses too low for them to 
evolve off the main sequence during the lifetime of the burst. The true gas 
usage timescale is thus 



The values of both /<ji s and / rem depend strongly on At. For example, for 
a burst duration At = 10 7 yr, only stars with M ^ 13 M@ have evolved off 
the main sequence. For a Miller-Scalo (1979) IMF with mi = 10, m u = 62, 
/rem ~ 0.4; if the lower mass limit is reduced to mi = 5, / rcm ~ 0.7. If 
e* = 0.1, / r em = 0.5, and f^ s = 0.5 (as is suggested by observations in the 
case of M82: see §4.1), then / = 5, and gas is effectively consumed at a 
much greater rate than that at which it is formed into stars. It must be 
emphasized again, however, that the effective values of / rem and f^ s are 
uncertain, depending not only on the duration of the burst but, for /d; s , on 
the physical conditions in the ISM, since this will presumably determine 
the recycling time for the gas. 



1 — e* 



ZdisM, + f rcm M s = fM s 



(11) 




(12) 
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3.2. AGN: 

The fueling requirements for AGN-powered ULIRGs are comparatively triv- 
ial: for the required rate M ~ 1 M yr _1 , the lifetime of the activity is 

— ^Mi^k- (13) 

where /i oss > 1 takes account of mass that is lost without accretion in the 
form of winds, etc. Of course, it is exceedingly unlikely that the gas reservoir 
is simply quiescently orbiting the black hole, waiting to be accreted. 

4. Superbubbles and Superwinds 

4.1. STARBURST: 

Massive stars input both energy and momentum into the ISM, through 
radiation, stellar winds and supernova blast waves. In systems with high 
rates of star formation, this input can have dramatic impact on the dy- 
namics of the ISM. In standard stellar wind theory (Weaver et al. 1977), 
the wind produces a hot shocked bubble in the ambient ISM; the swept-up 
ISM at the boundary will collapse to a thin shell provided the density is 
high enough for cooling to be efficient. 

In regions with many young, massive stars, a much larger scale ver- 
sion of such a wind bubble may result, as multiple stellar winds and es- 
pecially supernovae contribute to the development of a single large cavity. 
The evolution of such "superbubbles" can be largely understood simply 
by extending stellar wind bubble theory, scaling upward to the mass and 
energy input from a starburst (Chevalier & Clegg 1985; MacLow k, Mc- 
Cray 1988; Tomisaka & Ikeuchi 1988). Scaled to a canonical supernova rate 
^sn ~ 1 yr _ \ the mechanical energy and mass input rates are 

^SxM-^) (_*=_) ergs- (14) 

m sn ~ io (t^pt) m © y^ 1 ( 15 ) 

For a uniform ambient density n Q , the radius, expansion velocity, internal 
pressure, and thermal energy of the superbubble are (e.g., Heckman et al. 
1996) 

U~l.l(^) 1/6 t?/ 5 kpc (16) 
y = 200 f^V 7 V /5 kms- 1 (17) 
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?4~5x W 6 L% 5 5 n 3 J 5 tj 4/5 cm" 3 K (18) 

E th ~ 4 x 10 57 L 4 3.5*7 erg. (19) 
Finally, the soft X-ray (RO SAT-band) luminosity is approximately 

t lr .42 r 33 / 35 17/35.19/35 -l , on s 

L x ~ 10 L 43 ; 5 n D 7 t 7 erg s ; (20) 

this includes the effects of evaporation into the interior of the bubble. In 
these equations the mechanical luminosity of the OB association (approx- 
imated as constant) is L SN = 3 x 10 43 L43.5 erg s _1 and the age of the 
superbubble is t = 10 7 i7 yr. 

The expansion velocity (17) is generally supersonic, so the superbubble 
drives a shock into the surrounding ISM. The resulting shock luminosity 
depends on both the shock velocity and the ambient density: 

L sh « 2 x 10 40 Vf 00 n o A sh erg s 4 ~ 2 x 10 42 L 43 . 5 erg s^ 1 (21) 

where the shock velocity T4h = lOOVioo km s~ 4 and the shock area A 8 h is 
in kpc 2 (Dopita & Sutherland 1996); the second line assumes V s h is equal 
to the bubble expansion speed and uses equations (16) and (17). For shock 
speeds in the range V^h ~ 10 2 — 10 3 km s _1 , the resulting Ha luminosity 
L Ha ~ O.OlLshV^' 6 . The [O III] A4959 + A5007 line luminosity (including 
the contribution from the UV precursor) is comparable to the Ha luminos- 
ity at Vioo ~ 2 and gets relatively stronger as the shock speed increases, 
with the ratio scaling approximately as V^' 3 (Dopita & Sutherland 1996). 

The finite scale height of the gas disk means that the superbubble 
can expand entirely out of the disk, provided it is energetic enough. Such 
"blowout" will occur if the dimensionless parameter 

V = i^L ~ 10 4 ( -JL_Y 2 ( P ^ ) V „i/2 > 100 (22) 
L ISM VlOOpc/ ^10 7 cm- 3 Ky K ' 

(MacLow & McCray 1988) where L ISM = P ISM H 3 /t^ yn is the "luminosity" 
of the ISM over the dynamical timescale tdym defined as the time for the 
bubble to reach a radius of comparable to the scale height H of the gas layer 
(cf. equation [16]). Once the superbubble shell reaches an altitude Z £ a 
few H, it will begin to accelerate along the density gradient and fragment 
due to Rayleigh- Taylor instabilities. Provided this occurs on a timescale 
less than the lifetime of the OB association powering the supershell, the 
continuous injection of energy and momentum from the association will 
lead to the development of a galactic wind. 

There is extensive evidence for such large-scale galactic "superwinds" 
(e.g., Heckman, Armus Sz Miley 1990; Veilleux et al. 1994; Heckman et al. 
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1996), from morphologies, kinematics, and the presence of density /pressure 
profiles in reasonable agreement with the predictions of wind models (al- 
though there is not always agreement on the latter: see Veilleux et al. 1994). 
In the prototypical ULIRG Arp 220, the superbubble does not appear to 
have blown out, as the morphology and kinematics suggest a confined, 
shocked bubble. The starburst galaxy M82, on the other hand, appears to 
possess a freely expanding wind. Suchkov et al. (1996) argue that the wind 
in M82 must be mass-loaded, i.e., the mass flux in the wind is enhanced by 
a factor of ~ 5 over the SNe mass deposition rate, by entrainment or evap- 
oration of interstellar material. The implied mass lost through the wind is 
M\y ~ 5 x 10 7 tj Mq, which is comparable to the total molecular gas mass 
present in the nucleus. If this is typical of superwinds, then the fraction of 
available gas which is lost from the system is large, ~ 0.5. 

4.2. AGN: 

The possibility that AGN might drive powerful nuclear winds dates back 
at least to Weymann et al. (1982). More recently, this idea has been re- 
vived by Smith (1993, 1996) to explain the kinematics of the NLR and as 
a mechanism for cloud confinement. In a very provocative paper, Bicknell 
et al. (1998) have proposed that the radio jets in Seyfert galaxies are actu- 
ally dominated by thermal plasma, with mechanical luminosities of order 
-^mech ~ O.lLboi; the winds start off only mildly relativistic. Bicknell et al. 
make specific application to the excitation of the NLR in NGC 1068, and 
find additional support for this model from study of regions of jet-cloud 
interaction within the NLR. Interestingly, in a study of soft X-ray emission 
from large-scale nuclear outflows in a sample of nearby, edge-on Seyfert 
galaxies, Colbert et al. (1998) conclude that the outflows must be domi- 
nated by thermal plasma. This is an idea which is clearly deserving of far 
more study. If correct, it would replace the still unanswered question "Why 
are some objects radio-loud and others radio-quiet?" with another: "Why 
are some jets dominated by nonthermal plasma?" 

5. Global Energetics of the ISM 

Most ULIRGs are inferred to possess large quantities (M ^ 10 9 Mq) of 
high-density gas (n H ~ 10 4 cm -3 ), characterized by large velocity disper- 
sions (cr up to ~ 100 km s _1 ); see e.g., Tacconi, this volume. This implies 
that the energy dissipation rate within the ISM is very large: dimensionally, 




(23) 
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where R = 100i?ioo pc, <J = IOO07 km s 1 , and the coefficient 77 ^ 1 (e.g., 
MacLow 1998). Since the total non-rotational kinetic energy of the ISM is 

£ kin ~ \m<? ~ 10- [JL_) a * erg , (24) 

the energy loss timescale is 

. -^kin _ w 1n 5^l00 /nr\ 

idiss ~ -7; — ~5x 10° yr. (25) 

#kin W7 

Equivalently, since the inferred area filling factor of dense ISM in ULIRGs is 
close to unity (e.g., Scoville, Yun, & Bryant 1997), the energy loss timescale 
(assuming strongly dissipative collisions) tdiss ~ R/a ~ 10 6 -Rioo/ CJ 7 Y r ) 
essentially identical with the above estimate. Since the energy dissipation 
timescale is much shorter than the probable lifetime of the ULIRG, the 
random bulk motions of the ISM must be continuously powered. 

5.1. STARBURST: 

As noted in §4, the mechanical energy input from supernovae for a L = 
10 12 Li2 L Q ULIRG is 

^-^(r^GoSy (26) 

which is adequate to supply the necessary kinetic energy provided that a 
substantial fraction is absorbed by the dense ISM (and that momentum 
transfer to the clouds is reasonably efficient); given the large area covering 
factors inferred in ULIRGs (e.g., Scoville et al. 1997), this seems quite 
reasonable. 

5.2. AGN: 

Supporting the mechanical luminosity of the ISM requires either (1) rel- 
atively efficient conversion of radiation to mechanical energy, or (2), me- 
chanical (momentum) input from an AGN wind, as discussed above, with 

^mech ~ 0.1%-Lbol- 

6. Dynamical Support by Radiation Pressure 

As we have just seen, there is enough energy available from either a star- 
burst or an AGN to maintain the energy loss rates estimated for the ISM 
in ULIRGs. A related topic is the momentum input rate. Scoville et al. 
(1995) suggested that radiation pressure might contribute significantly to 
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the support of the gas against gravity in the z = 2.3 ULIRG F10214+4724, 
in order to reconcile the CO-derived gas mass with a dynamical mass (es- 
timated from the CO source size and line width) which is an order of 
magnitude smaller. 

Assume that a cloud at radius r absorbs all of the incident flux (and 
therefore photon momentum) incident upon it from the central source of 
radiation. In equilibrium, the condition that the radiation provide enough 
pressure to support the cloud against gravity is just 

GMnm H N H = L 
r 2 Airr 2 c 

where M is the total mass enclosed within radius r, fi is the mean atomic 
weight in units of the hydrogen mass m H , and N H is the hydrogen column 
density through the cloud. Simple rearrangement shows that the cloud col- 
umn density which can be supported is 

JV2*~— i: ~ 3.5 x 10 22 (^M.) cm -2 (28) 

which depends on the luminosity to mass ratio, here expressed in solar 
units. (Note that although I have done the calculation by balancing the 
radial gravitational attraction against the radiation pressure from a central 
source, the identical result is obtained by calculating the column which 
can be vertically supported above a radiating slab, which may be a more 
relevant geometry for starburst-powered ULIRGs). 

The limiting column density which can be supported is much smaller 
(by a factor 10 — 100) than the area-averaged molecular column densities 
which have been derived in the nuclei of ULIRGs, for which the values of 
L/M generally lie in the range 10 2 — 10 3 , suggesting that radiation pres- 
sure support is not important. The expression for the critical column (28) 
assumes that the photons give up all of their momentum immediately, in 
other words, it has been assumed that the ISM is optically thin to the 
re-radiated photons (generally in the mid- to far- infrared) . If the ISM is 
optically thick to these photons, then multiple scattering of the radiation 
will occur. However, this will not increase the column which can be sup- 
ported. Assume that the medium is very optically thick, so that a large 
number of scatterings occur. In this case the radiation field will become 
nearly isotropic. The degree of anisotropy (which determines how much 
work the radiation field can do against gravity) is then set by the net flux 
outward; this is (in spherical geometry) precisely L/AirR 2 , and so the value 
of N™ d is unaltered. 
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7. Pressure in the ISM 

In any model for ULIRGs, an inevitable result is that the pressure in the 
ISM must be extremely high: 

• Simply from the derived gas densities and observed turbulent veloci- 
ties, the inferred turbulent pressures must be high: 

5x108 (i^)G^) 2 -- 3k - < 29 > 

• If the nuclear (R ~ 100 pc) ISM essentially radiates like a blackbody, 
as proposed by Downes & Solomon (1998), then the radiation pressure 
in this region 

p - = i~i^k < 3 »> 

(where u is the radiation energy density) and so therefore 

P rad - 3 x 10 8 L 12j Rr 2 cm" 3 K . (31) 

• In an AGN-powered ULIRG, the pressure of radiation from the central 
source of luminosity will force the gas pressure to large values: 

P rad ~ 8 x 10 7 (J^j L 12 R W 2 cm- 3 K (32) 

where / a b s is the fraction of the incident radiation which is absorbed 
in a thickness 5r <C r c , the cloud radius (see Barvainis et al. 1997). 

• In a starburst-powered ULIRG, with enough energy input from SNe 
to drive a superwind, the pressure in the injection region (i.e., on the 
spatial scale of the starburst) is 

P ~ O.lM^Ll^R-i (33) 

(Chevalier & Clegg 1985) where R SB is the radius of the starburst 
region; hence 

- 3 x 108 (w) ( lo^Fi) ( Wi;) " ■ < 34 » 

All of these estimates agree to order of magnitude, implying that the 
pressure in the ISM P ISM ~ 10 8 cm -3 K in ULIRGs. This is 3 — 4 orders of 
magnitude larger than the pressure in the average ISM in the solar neighbor- 
hood; only the densest cloud cores approach such high pressures. Similarly 
high pressures have been inferred in more modest starburst galaxies such 
as NGC 1808 and NGC 3256 (e.g., Aalto et al. 1994, 1995). 
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8. Local Energetics: Energy Input to Molecular Clouds 

8.1. STARBURST: 

The power input from an OB association in the form of radiation, stellar 
winds, and supernovae is 

£ OB ~ 6 x 10 35 iV* erg s" 1 (35) 

(McCray & Kafatos 1987) where iV* is the number of stars in the association 
with M > 8 Mq. For lower and upper stellar mass cutoffs of mi = 1 and 
m u = 60 Mq, respectively, (which gives 1 SNe for every ~ 100 Mq turned 
into stars) the energy input to the cloud from which the association formed 
over the cloud's lifetime At is 

where M c is the cloud mass and e* is the star formation efficiency as before. 
The cloud binding energy (for mean cloud density n H = 10 4 ri4 cm -3 ) 

the ratio of energy input to binding energy is thus 

Eob nK ( U\( M c \- 2 / 3 _i/3 ( At 



E B - 25 {m){wwJ "'{my,) < 38 » 

which depends relatively weakly on the assumed cloud parameters. Thus, 
in principle, only a few percent of the available power needs to be input 
to the cloud in order to disperse it (see e.g., Williams & McKee 1997). 
However, most of this energy will of course be radiated away. 

8.2. AGN: 

If the clouds in the ISM have a direct view of the active nucleus, then the 
integrated energy input (assuming At < the lifetime of the AGN) is simply 



Eagn ~ ^W AclM 



where A c \ is the projected area of the cloud. Thus the ratio of AGN energy 
input to cloud binding energy is 

E AGN in6r p _ 2 , ( M c n- 1 _i/ At \ ,, n . 
- 10 L 12 iW abs —c— n 4 (40) 



E B 1Z W0JaDS \ltf MqJ 4 Vl0 7 yr 



12 



PHILIP R. MALONEY 



which is several orders of magnitude larger than the energy input calculated 
above for a starburst. (Since the bolometric luminosities are the same in 
both cases, the same amount of radiation - mostly re-radiated in the far- 
IR - must of course be incident on a typical cloud in a starburst-powered 
ULIRG. Precisely because this radiation is dominantly at infrared wave- 
lengths, however, it will have far less impact on the ISM in a ULIRG than 
the radiation from an AGN, serving mostly to increase the typical ISM 
temperatures.) Evaporation and ablation of clouds which are exposed to 
the AGN can be very significant, leading to destruction of the molecular 
ISM (Begelman 1985). 

Because of the large gas column densities inferred in the nuclei of 
ULIRGs, however, it is likely that most of the ISM does not have an un- 
obscured view of the AGN. Here the fact that the spectrum of an AGN 
extends to high energies {E ~ 100 keV) is of great importance, as these 
photons can penetrate large column densities before being absorbed. This 
means that the amount of energy deposited by X-rays can be very large 
even if there is a significant column between a cloud and the radiation 
source: the energy deposition rate per unit volume is approximately 

nH x ~ 3 x 10- 17 (A) L 12 R W 2 ( 102 f a c ^_ 2 ) 1 n 4 ergcm- 3 s- 1 (41) 

(Maloney, Hollenbach, & Tielens 1996, hereafter MHT), where f x is the 
fraction of the bolometric luminosity emitted in the 1 — 100 keV energy 
band, and N at t is the column density of neutral hydrogen between the cloud 
and the AGN. For shielding columns N at t ~ 10 22 , only hard X-ray (E > 
1 keV) photons will be present. Unlike the case of the photodissociation 
regions (PDRs) powered by stellar UV photons, where dust absorption 
leads to exponential attenuation of the heating and ionization rates with 
depth, the XDRs powered by AGN are characterized by an only linear 
decrease of the energy input with increasing column density (the exact 
scaling with column density depends weakly on the spectral index of the 
incident radiation field: MHT). The total energy input to the cloud is then 

~ 1055 (It) L " « iwk) (wvr) «* < 42 » 

where the shielding column has been scaled to A^ att = 10 22 A?22 cm -2 . Com- 
pared to the cloud binding energy, 

^agn ,rx(fx\ T p-2 v -i / 1O 5 M \ 2/3 _ 1/3 / At \ 

-gr- ~ 10 [ U ) L 12 R 100 N 22 j n 4 (^-j erg. (43) 

Thus, even if the shielding colummn between the X-ray source and a cloud 
is iV a tt ~ 10 24 cm -2 , the X-rays from the AGN can still dump ~ 10 3 times 
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the cloud binding energy into a cloud over its lifetime. As in the case of 
a starburst-powered ULIRG, most of this energy is simply radiated away 
(although photoevaporation will occur if the clouds are heated enough to 
raise the sound speed above the escape speed). It does mean, however, that 
it is extremely easy to "light up" the interstellar medium in a ULIRG, 
simply in consequence of the energy input rate, whether it is powered by 
an AGN or a starburst. 



9. X-Ray Irradiation in AGN-Powered ULIRGs 

As alluded to above, the large column densities which hard X-ray photons 
can traverse before being absorbed, combined with the large energy per pho- 
ton, allow X-ray luminous AGN to have a profound impact on the physical 
and chemical state of the ISM in AGN-powered ULIRGs. In this section I 
sketch briefly some of the results and their observational consequences. 



9.1. PHYSICAL AND CHEMICAL STRUCTURE 

Detailed models of the physics and chemistry of XDRs have been discussed 
by MHT. The important parameter is H x , the X-ray energy deposition 
rate per hydrogen nucleus. The heating rate (as noted earlier) and the 
molecule destruction rates scale as nH x , whereas the cooling and molecule 
formation rates are generally oc n 2 ; thus the physical and chemical state 
scales approximately with H x /n. It is convenient to express this in terms 
of an X-ray ionization parameter (the ratio of photon flux to gas density) , 

L x 4irF x 

« = ^ = — < 44 > 

where F x is the X-ray flux incident on the face of the cloud. For opti- 
cally thin gas, the ionization parameter determines the physical state {e.g., 
Tarter, Tucker & Salpeter 1969). In the case of interest here, where the 
clouds are optically thick to the X-rays (up to some maximum energy), the 
physical conditions depend on an effective ionization parameter, 

AttF x L aa 

( " = ^°- 1 ^KM (45) 

where the 1 — 100 keV X-ray luminosity is L x = 10 44 L44 erg s~ x and, as 
noted earlier, the scaling index a ~ 1 for typical AGN spectral indices. 
This effective ionization parameter is directly proportional to H x /n, with 
H x /n ~4x 10~ 25 £ cff erg cm 3 s _1 . 

Figure 1 shows the physical and chemical structure of an XDR as a 
function of £ e ff or, equivalently, distance from the X-ray source. This cal- 
culation assumes L x = 4 x 10 44 erg s _1 , or f x « 0.1 for Lboi = 10 12 L®. 
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Figure 1. XDR physical conditions as a function of radius (or £ e ff), for L x — 4 x 10 44 
erg s~\ iVatt = 10 22 cm~ 2 , and n = 10 5 cm" 3 . 

The gas density is constant at n = 10 5 cm -3 , and the column between the 
X-ray source and the gas is iV att = 10 22 cm~ 2 . (Cloud column densities 
N c \ = 2 x 10 22 cm -2 and linewidths Ay = 5 km s _1 have been assumed in 
calculating the cooling of optically thick species.) Close to the X-ray source, 
the gas is warm, T ~ 10 4 K, atomic, and weakly ionized, with ionization 
fraction x e ~ 0.01 — 0.1. With increasing distance (or increasing column 
to the X-ray source) the temperature and ionization fraction decreases and 
the abundances of molecular species increase; at the lowest values of £ e ff> 
the temperature has dropped to T 25 K, and nearly all of the gas-phase 
carbon is in CO. Atomic oxygen is abundant over the entire range of £ e ff- 

9.2. MOLECULAR SURVIVAL 

Figure 1 shows that, for a gas density n = 10 5 cm -3 , CO only becomes an 
abundant species for R £ 100 pc, whereas many ULIRGs show bright CO 
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emission on smaller spatial scales. What densities are needed for gas to be 
molecular in an AGN-powered ULIRG? 

Define the fractional abundance of a species by Xj = n^/n, where rii is 
the number density of species i and n is the total density of hydrogen nuclei, 
n = n H + 2n H2 ; note that this definition means that the molecular hydrogen 
fraction x H2 = 0.5 if all of the hydrogen is in molecular form. In table 1, I 
give the density n m required in order for the CO abundance Xco = 

10" 4 for 

three different distances, R = 10, 100, and 1000 pc from the nucleus, and 
three different values of the shielding column, iV att = 10 22 , 10 23 , and 10 24 
cm -2 . Also listed in the table are the corresponding molecular hydrogen 
fraction x Ha and the gas temperature T. 

TABLE 1. Densities for CO Survival 
iVatt n m x H2 T(K) 



R = 10 pc 



10 22 


1.5 


X 


10 8 


0.07 


903 


10 23 


2.5 


X 


10 7 


0.10 


694 


10 24 


4.8 


X 


10 6 


0.13 


545 








R = 100 pc 






10 22 


1.9 


X 


10 6 


0.14 


409 


10 23 


3.3 


X 


10 5 


0.19 


233 


10 24 


8.5 


X 


10 4 


0.27 


146 








R = 1000 pc 






10 22 


3.7 


X 


10 4 


0.30 


125 


10 23 


4.5 


X 


10 4 


0.47 


74 


10 24 


1.4 


X 


10 4 


0.48 


42 



Table 1 clearly demonstrates how powerful the effect of an X-ray lumi- 
nous AGN on the ISM can be. At a distance of 100 pc from the nucleus, 
even with a column of iV at t = 10 24 cm~ 2 in front of the X-ray source, the 
gas density must be few x 10 4 cm 3 in order for CO to be abundant: 

for a density n = 2 x 10 4 , the CO abundance x co < 10~ 6 . The temperatures 
are considerably higher than the typical ISM in the Milky Way; note also 
that the thermal balance has been calculated including only the heating 
due to X-rays. Clearly, high gas densities are mandatory for molecular sur- 
vival in AGN-powered ULIRGs, even when the AGN is buried beneath a 
large column density of gas and dust. 
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10. Far-Infrared Fine-Structure Lines 

In the post-ISO era, surely no one needs to be reminded of the wealth 
of information provided by far-infrared fine-structure line observations, a 
point that was made repeatedly throughout this meeting. As discussed by 
MHT, XDRs are copious sources of emission in these lines. 

Figure 2 shows the emergent surface brightnesses in a number of impor- 
tant mid- and far- infrared lines, for the same XDR model parameters as 
used for Figure 1. Also plotted is the total far-infrared continuum surface 
brightness from grain emission. Of particular note are: (1) the extremely 
high line surface brightnesses which can be produced; (2) major cooling 
lines such as [O I] 63/jm and [Si II] 35//m are bright over nearly the entire 
parameter space; and (3) compared with the far-IR continuum produced 
in the XDR, the surface brightnesses in these same lines can be extraor- 
dinarily high, with line-to-continuum ratios approaching 10%. In contrast, 
the line-to-continuum ratios produced in PDRs are usually of order 0.1%, 
and rarely exceed 1% (e.g., Tielens &: Hollenbach 1985). This difference is 
a reflection of the fact that in PDRs the gas heating is limited by the grain 
photoelectric heating efficiency; in contrast, in XDRs about as much energy 
is deposited into the gas as into the grains, so that major cooling transi- 
tions can carry a large fraction of the total deposited energy (see MHT for 
details). The other important point to note is that all of these lines arise 
from neutral or singly-ionized species; because of the large column densities 
which X-ray photons can traverse, the volume of the XDR can be vastly 
larger than the volume of any high-ionization region produced by the AGN. 

10.1. PARTICULARLY INTERESTING LINES 

A number of other infrared transitions are of particular interest for ULIRGs. 
The near-infrared vibration-rotation lines of molecular hydrogen have been 
detected from a number of ULIRGs (most spectacularly in the not-quite- 
an-ULIRG NGC 6240); the 1.64/xm fine-structure line of [Fell] is frequently 
detected as well. Observations with ISO have produced detections of the 
[Nell] 12.8/um line in a number of ULIRGs as well, and the first extragalactic 
detections of the pure rotational transitions of H2 , including observations 
of the archetypical ULIRG Arp 220 (Sturm et al. 1996). 

In Figure 3 I have plotted the emergent surface brightnesses for the 
v = 1 — and v = 2 — 1 S(l) lines of molecular hydrogen, the 1.64/im [Fell] 
line, Br7, and the [Ne II] 12.8^m line. As in Figure 2, the high line surface 
brightnesses are striking. The extent of the [Fe II] emission is limited by 
gas temperature, since the upper state lies E/k ~ 10 4 K above the ground 
state. The drop-off in molecular hydrogen emission at log£ c g- « 2 is also 
due to declining T (the v = 1 and v = 2 levels lie approximately 6000 and 
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Figure 2. Emergent surface brightnesses in important mid- and far-infrared lines for 
the XDR model of Figure 1. 

12,000 K above ground, respectively); the secondary rise is due to direct 
nonthermal excitations of the levels. For a fixed density, the [Fell] emission 
always lies interior to the molecular hydrogen emission (see also Maloney 
1997). The run of density with £ c g- (or radius) will determine the ratios of 
the H2 and [Fe II] lines with respect to each other and to Br7. 

An additional interesting result is the intensity of the [Ne II] 12.8/um 
line. This is the single most important coolant over a substantial fraction 
of the plotted parameter space; as much as ~ 15% of the deposited X- 
ray energy emerges in this line alone. The resulting surface brightnesses 
are extraordinarily high, and the total line luminosity can be a substantial 
fraction of the total hard X-ray luminosity. 

The pure rotation lines of molecular hydrogen are of particular interest 
because they are much less subject to extinction than the near-infrared 
lines, and because they provide a tracer of warm (few hundred K) molecu- 
lar gas, unlike the vibration-rotation bands which require T ^ 2000 K for 
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Figure 3. Emergent surface brightnesses in several important near- and mid-infrared 
lines for the XDR model of Figure 1. 

effective excitation. The S(0), S(l), S(3), and S(5) lines have been detected 
in NGC 1068 (Lutz et al. 1997) and the S(l) and S(5) lines (at 17 and 
6.9/t/m, respectively) have been detected in Arp 220 (Sturm et al. 1996). 

In Figure 4 I have plotted the surface brightnesses of the S(0), S(l), S(3), 
and S(5) lines for the XDR model of Figure 1. Current observations lack 
the spatial resolution to do more than measure the total line luminosities. 
However, even this should prove extremely useful for constraining models 
of the ISM in ULIRGs, and XDR models in particular (Lutz et al. 1997). 

10.2. THE HCN/CO RATIO 

In ULIRGs, the ratio of the HCN J = 1 -> line to that of CO is typically ~ 
0.2 (Solomon, Downes, & Radford 1992), an order of magnitude larger than 
the value typically observed in normal spirals. In the nuclear (R ;$ 100 pc) 
region of NGC 1068, this ratio reaches unity (Tacconi et al. 1994). Such high 



o 
i- 



tafi 

o 



THE ISM IN ULIRGS 
Log R (pc) 
1.5 2 2.5 3 3.5 

| — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 



19 



-2 - 



-4 - 



-6 



n = 10 5 cm -3 




-8 - S (!) 



_i i i i_ 



_i i i i_ 



-10 

-1 -2 -3 -4 -5 

L °g £eff 

Figure 4- Emergent surface brightnesses in the S(0), S(l), S(3) and S(5) pure rotation 
lines of Hb, for the XDR model of Figure 1. 



values have been interpreted as meaning that the density is high (n ^ 10 5 
~ 3 ), so that both the HCN and CO lines are approximately thermalized, 



cm 



and that both of the transitions are optically thick, or nearly so. However, 
the inferred high densities must be regarded with some caution, as the 
HCN rotational transitions can be pumped through a bending mode via 
absorption at 14/^m; this process is plausibly important in both starbursts 
and AGN (Aalto et al. 1995; Barvainis et al. 1997). 

The line center optical depth for a transition at frequency v is 



t„ = 



8vrz/ 3 



A 



9 j Ni 



1 



(46) 



where Aji is the Einstein A— coefficient, j and i denote the upper and lower 
states, the gs are the statistical weights, Ni is the column density in the 
lower state, AV is the linewidth, and T cx is the excitation temperature 
characterizing the relative populations of levels i and j. Assuming equal 
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excitation temperatures for the two transitions and that hv/kT ex <C 1, the 
ratio of r D = 1 columns for the J = 1 — > transitions is then 

« = * 1.8 x 10-. (47) 

iV (CO) |/X H cn| 2 ^hcn 

This ratio is much less than unity because the ^4— coefficient for the HCN 
transition is so much larger than for the corresponding transition of CO, 
in consequence of its much larger dipole moment ji. Thus, if the HCN 
abundance is more than a fraction of a percent of the CO abundance, the 
HCN J = 1 — > transition will be optically thick if the CO transition 
is. Of course, if the CO J = 1 — > optical depth is very large, then a 
much smaller HCN abundance will still lead to an optically thick HCN 
line. Furthermore, it is very important to note that equation (47) gives 
the ratio of required column densities in J = 0. Because of its much lower 
critical density, it is likely that the partition function for CO will be close 
to the LTE limit for the relevant conditions, while that of HCN will not, 
so that the CO molecules will be distributed over a much larger number 
of rotational states. Hence the ratio N /N tot will be much larger for HCN 
than for CO, also reducing the required HCN abundance. 

In the quiescent interstellar medium in the solar neighborhood, the op- 
tical thickness condition is not fulfilled, as the HCN abundance is usually 
^hcn ~ lCr 4 x co - Furthermore, there are indications that in very active 
galactic nuclei (either high SFR or AGN), the optical depths in the CO 
lines are only modest (Aalto et al. 1995; Barvainis et al. 1997), so that 
the HCN lines are not optically thick in reflection of very optically thick 
CO lines. However, the increased ionization rate in an XDR can lead to an 
enhanced abundance of HCN relative to CO (Lepp & Dalgarno 1996). In 
Figure 5 I show the abundances of CO, HCN, and the ratio x UC n/x C o as 
a function of £ e ff and distance. The AGN parameters are the same as pre- 
viously, but I have increased the density to n = 10 8 cm -3 simply to make 
the gas dominantly molecular within 100 pc of the nucleus. The HCN/CO 
ratio is greater than 10 -3 over about a factor of three in radius. Thus in 
an XDR the HCN abundance can be elevated sufficiently relative to CO to 
explain the observed line ratios. (In hot star-forming cores, the HCN/CO 
ratio also reaches values of order 10 -3 [Blake et al. 1987], so elevated values 
may be expected in starburst-powered ULIRGs as well.) 



11. PAH Emission 

The mid-IR emission features attributed to polycyclic aromatic hydrocar- 
bons (PAHs) have been used to discriminate between AGN and starbursts 
in ULIRGs (e.g., Lutz et al. 1998), based on the observation that these fea- 
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Figure 5. HCN and CO abundances and their ratio as a function of £ e ff- The assumed 
AGN and cloud parameters are the same as Figure 1, except that the assumed density 
has been increased to n = 10 s cm~ 3 . 



tures are weak or absent in classical AGN but generally strong in starbursts. 
Does the presence of PAH features imply the absence of an AGN? 

Voit (1992) has investigated the effects of X-ray irradiation on PAHs. 
Absorption of an energetic photon leads to destruction of the PAH, through 
two processes: (1) photo-thermal dissociation {i.e., evaporation), and (2) 
Coulomb "explosion", in which the PAH is doubly ionized by the incident 
photon and subsequently fragments as the repulsive Coulomb force effec- 
tively reduces the binding energy of the PAH. To estimate the magnitude of 
the destruction timescale, assume that every absorption of an X-ray photon 
by a PAH leads to its destruction 1 : 

Txd ~ 2 (ol) Li2jR ioo^22 6 yr (48) 

x This is likely to overestimate the PAH destruction rate; far more laboratory work 
{e.g., Jochims et al. 1996) needs to be done in this area. 
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assuming a flat (yL v = constant) spectrum. Even if iV a tt = 10 24 cm -2 , the 
PAH X-ray destruction time is only r x d ~ 3000 years at 100 pc. 

To determine whether PAHs can survive in the face of such short de- 
struction timescales, we need to compare r x( j with the rate at which the 
PAHs re-accrete material (this is dominated by accretion of C + if the PAHs 
are neutral). Since a photo-destruction rate is balanced against a density- 
dependent reaction rate, this again leads to a critical ionization parameter 
at which r x d = r acc , which, using the results of Voit (1992), is 

^M. 5xl0 -»(^)</ PAH ^ ,49) 

where Y is an effective sticking coefficient for accretion, the number of 
carbon atoms in the PAH is = 50A%o, and / PA h° is the fraction of 
PAHs which are neutral. S(A r H ) is a dimensionless, spectrum-weighted cross 
section for photodissociation of carbon from a PAH, (i.e., the photon flux- 
weighted photodissociation cross-section, normalized to the PAH geometric 
cross-section), approximately given by ^(N H ) ~4x 10~ 5 /A^22 (see Figure 3 
in Voit 1992). This gives the critical ionization parameter for PAH survival 
as 

U cr ~ 1.3 x 10~ 4 (j^J NH 2 f PAll of c+ N 22 (50) 

where f c + = x c +/x c (total). From this expression, we find that PAHs will 
survive against X-ray destruction at a radius (with f c + = /pah° = 1) 



R cr ~ 100 



0.1 J VlOOK 



1/2 / p \ -V2 



Note the importance of large (N22 ^ 1) shielding columns and high pres- 
sures for the survival of PAHs. Since, as we have seen, both large column 
densities and high pressures are characteristic of the ISM in ULIRGs, it is 
quite probable that PAHs will in general resist destruction by X-ray irradi- 
ation even in AGN-powered ULIRGs, and so the presence of PAH features 
should not be taken as conclusive proof of the absence of an AGN. 
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